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Abstract 

The perfect fluid in the context of a covariant variable speed of light theory 
proposed by J. Magueijo is studied. On the one hand the modified first law of 
thermodynamics together with a recipe to obtain equations of state are obtained. 
On the other hand the Newtonian limit is performed to obtain the nonrelativistic 
hydrostatic equilibrium equation for the theory. The results obtained are used to 
determine the time variation of the radius of Mercury induced by the variability 
of the speed of light (c), and the scalar contribution to the luminosity of white 
dwarfs. Using a bound for the change of that radius and combining it with an 
upper limit for the variation of the fine structure constant, a bound on the time 
variation of c is set. An independent bound is obtained from luminosity estimates 
for Stein 2015B. 



1 Introduction 

There are several very different motivations for studying the variation of fundamental 
constants. The coincidence of large dimensionless numbers arising from the combina- 
tion of different physical constants led Dirac to propose the large number hypothesis 
and predict a time variation of them [1] , [2] ■ Theories with varying constants can also 

*racker@ecni. ub.es - Present address: Departament d'Estructura i Constituents de la Materia and 
ICCUB, Institut de Ciencies del Cosmos, Universitat de Barcelona, Diagonal 647, 08028 Barcelona, 
Spain. 

^ sisterna@mdp . edu . ar 

^ vucet ichQfcaglp . unlp . edu. ar 



1 



be a way of implementing Mach's principle. Extra-dimensional theories like superstring 
or Kaluza-Klein theories reduce in the low energy limit to effective theories in which 
the fundamental constants may vary in space and time. 

Although no variation has been found in most experiments and observations per- 
formed up to date, the results from analysis of spectra from high redshift quasar ab- 
sorption systems remain controversial. Some works have reported a variation in the 
fine structure constant [3], [1], [3], [B], [T|, [S], while other studies give null results [H], 
[To] . [TT] . [T2] . Besides the motivations mentioned above, variable speed of light (VSL) 
theories are interesting because they could solve several cosmological puzzles [13], [H] . 

m- 

In this work we study the thermodynamics and Newtonian limit of the varying 
speed of light theory developed by J. Magueijo [H]. This theory was an improvement 
of a previous version p3], in order to account more carefully of local Lorentz invariance 
and the dynamics of the scalar field behind c's variation. The theory here considered 
has been strongly criticized in |T7], responded in [18j, criticized again in [19], [20], and 
answered back in [21]. As we point out at the end of the paper, in Magueijo's theory 
the interaction between matter and the scalar field is not explicitly shown in full detail. 
This can be seen in the effective quantum creation of particles (see Sec. VI of Ref. |16]). 
even in the absence of an explicit interaction term in the matter Lagrangian (6 = case; 
see below). Nevertheless, being a theory that gathered a lot of attention in the recent 
past, we consider that it is very suitable as a first application of the general framework 
developed in this work, which we plan to apply to other (formally less controversial) 
theories such as bimetric theories, in future publications. 

After a summary of this VSL theory (Sec. 2) and of a Lagrangian approach to de- 
scribe perfect fluids (Sec. 3), the first law of thermodynamics and a recipe for obtaining 
equations of state are derived (Sec. 4). It is shown that the field associated with the 
variation of the speed of light can formally be considered as a new thermodynamic 
variable. Regarding this point we note that this field has two properties that are not 
common in thermodynamic variables: local universality and long scale variation. This 
has been remarked in a work in the context of another theory with variation of physical 
constants [22]. In Sec. 5 we perform the Newtonian limit and in Sec. 6 we apply all 
the results to study the evolution of the radius of Mercury and derive a bound on the 
time variation of c. Section 7 devotes to the estimation of the scalar contribution to 
the luminosity of a white dwarf, obtaining a stringent upper bound for the variation 
of c. In Sec. 8 we state our conclusions and we leave for the Appendix some results 
concerning the space and time dependence of the scalar field. 

2 Brief description of the VSL theory 

In the covariant and locally Lorentz invariant VSL theory proposed by Magueijo c plays 
three different roles: 

a) c is a dynamical field: 
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The spacetime variations of c are represented by an adimensional scalar field 
so that 

c = coe^ (1) 

where cq is a constant. 

The general relativity (GR) action is modified and becomes 

/ = j d^x (^e'^'f'iR -2A- ^V^ijWij) + ^^e'^£^^ . (2) 

The metric is taken to be 77^,^ = Diag{—1, 1, 1, 1). Cm is the matter Lagrangian 
and K, a, and b are three parameters of the theory. We will take a — 6 = 4 as 
in [16]. The matter Lagrangian is required to have no explicit dependence on c 
{minimal coupling condition) , which leads to the second role played by c. 

b) c parametrizes the variations of the other "constants" : The minimal coupling 
condition fixes the scaling with c of all the Lagrangian parameters up to the h{c) 
dependence, which is taken to be ^ oc c^~^, where q is the fourth parameter of the 
theory. For example, since the Compton wavelength appears in the Lagrangian 
of a quantum particle, the dependence on c of the mass will he m (x h/c oc c''"^. 
In a similar fashion it can be determined that the charge of a quantum particle 
(e), the Bohr radius (r;,), and the fine structure constant (a) are proportional to 
c'^, c~'^, and c'^, respectively. 

c) c is an integrating and conversion factor: The theory is covariant and locally 
Lorentz invariant in a generalized way explained in [16j. The key point is the 
use of an x° coordinate, instead of time, in all geometrical formulas. The main 
difference with constant c theories is that local measurements of space and time 
[dx and dt) are not generally differentials of a coordinate system (the partial 
derivatives do not commute). Nevertheless it is always possible to find integrat- 
ing factors such that dtip^ and dxtp^'^ are perfect differentials, with j3 another 
parameter which, however, will not appear in the equations of our work, because 
the calculations will be done using the x° coordinate. 



Varying the action with respect to the metric and ip leads to the equations: 



Ui, + aV,i,V^i. = jip^!^ (aT - 2ap,c^ - 26£„.) + i| , (4) 

where T'^'^ is the matter stress energy tensor and T is its trace. Besides, pA is the 
mass density corresponding to the cosmo logical constant, pA = and A is a linear 
combination of the matter and geometric cosmological constants, A = A + ^^A^ 
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After applying the Bianchi identities to Eq. and using Eq. (jl]), an equation for 
the divergence of T^'^ is obtained: 



= -^;.T%6 + ij.^bC^ - A^.^ . (5) 

Note that matter energy is conserved only when 6 = 0, in all other cases there is 
exchange of energy between matter and the field. 

The presence of the ip field can modify the law of conservation of the number of 
particles and the normalization condition for the four-velocity, e.g. it is found that for 
a classical particle f/^ = Uq (c/cq)"^ 7^ —1 [16j. Besides, the energy density and the 
total energy of a body in hydrostatic equilibrium will also vary if c is not constant (we 
must distinguish between energy density and total energy because the size of a body 
can change in time as a result of the variation of c). In addition, the c dependence 
of the mass is different for a classical and a quantum particle. Finally, the matter 
Lagrangian, which is not unique, appears in the equations of the theory. To take 
into account all these effects and ambiguities it is convenient to introduce four new 
parameters, qi,q2,q3, and q^. 

^g^^^^ _ ^91^2 _ ^ Generalized normalization 

1^ ^Yie four-velocity. ^ 

/ n^TT,y\ Generalized conservation 

inC^^U )■;. = r 1 1 ( ' ) 

' of particle number. ^ ' 

_ c dependence of the energy density of , . 

^ a body in hydrostatic equilibrium. ^ ' 

U — U ^ dependence of the total energy of 

a body in hydrostatic equilibrium. ^ ^ 

The " " subscript denotes the value of these quantities when = 0. 



3 The perfect fluid Lagrangian 

The task of obtaining a Lagrangian for the perfect fluid is not a trivial one due to the 
constraints imposed by the normalization of the velocity and the conservation of the 
number of particles. A. H. Taub gave one in 1954 [23] and another one was given by B. 
F. Schutz in 1970 [21], using a different but equivalent approach. Note that equivalent 
Lagrangians in usual theories (i.e. differing in a divergence) may not be equivalent 
in Magueijo's one, because of their explicit appearance in the evolution equations. In 
our work we have used Schutz's Lagrangian, hence we present in this section a brief 
summary of his approach and then we show how it can be used in the VSL theory 
under study. 
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3.1 Schutz's Lagrangian 

Schutz uses a formulation of relativistic hydrodynamics based on the utihzation of 6 
potentials to represent the velocity: 

= fj,-' (0,, + CP,. + Os,,) , (10) 

where /i and s are the specific enthalpy (enthalpy per unit mass) and the specific 
entropy, respectively. The physical meaning of the other potentials is also explored by 
Schutz. 

The perfect fluid action is 

«+^) (--lY'^c^^, 

where p is the pressure of the fluid. Then the perfect fluid Lagrangian is Cm = P- The 
following steps must be followed to vary the action: 

1. Choose an equation of state for the fluid and write it in terms of /i and s: p = 
p(/i, s). 

2. When varying the action make use of the thermodynamic relation dp = Pmdp — 
PmTds (where is the rest mass density). 

3. Define the four-velocity vector field U,y in terms of 6 scalar potentials: Uu — 

4. The normalization of U is taken into account before starting the variations. This 
is done expressing in terms of 0, /?, 6', s, and g^": /j,"^ = —g^" {4>^^ + 03^^ + Os^n) 
{4>,v + il3,u + Os^u)- So the independent variables are (j),C,P,9,s, and g^'^. Any 
quantity appearing in the action must be considered a function of these variables. 

The Euler- Lagrange equations become 

Gixu - [{p + p)U^Uy + pg^y] = , 

{PmUn;. = 0, 

C/% = 0, U'^e^^^T, U^'p^^^O, and U^^^^O. 
The stress energy tensor is 

rpni^ ^ 2 ^ (v^£^) ^ p^u^'u" + pg^" . (11) 

It is important to note that (p+p), U^,p, Pm, and T were defined in terms of 0, ^, /3, 9, s, 
and g'^'^ after performing the variations, using the equations 

p — p{n, s) (equation of state), 

dp\ rp^Zlf^ 

dp) Pm \ds 

(P + P)^ PmP: = {.(t>,v + iP,v + OS^^) 
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3.2 Use of Schutz's Lagrangian in Magueijo's theory 

Using Schutz's Lagrangian in the action ([2]) and varying it with respect to (f),^,P,9, 
and s we obtain 





(12) 


urs,, = 0, 


(13) 


we^, = T, 


(14) 




(15) 


U'iu = 0. 


(16) 



Moreover, varying Cm with respect to g^'^ leads to 

T^"" = {p + p)U^'U'' + pg^"" , (17) 

where p,p,Pm, and are defined from 0, ^,/3,6', and s exactly in the same way as 
in Schutz's theory. In the VSL theory they can depend on ip, but they will coincide 
with the usual quantities when ip = 0. Besides, with these definitions of p,p, and U^, 
the stress energy tensor in the VSL theory has the same form as the usual one for a 
perfect fiuid with energy density p and pressure p (both quantities being measured in 
a momentarily comoving reference frame). These facts make it reasonable to consider 
p,p,Pm, and f/'^ as the energy density, pressure, rest mass density, and four-velocity of 
the perfect fiuid in the VSL theory that is being studied. 

Note that by definition U^U^ = —1, so qi = 0. This is different from Magueijo's 
result for a classical particle, e^^W^U^ = constant, where the definition of the velocity 
is = ^ = Although there is no contradiction with this, one has to be careful 
with the interpretations given to U^. 

4 Thermodynamics 
4.1 First law 

When energy is conserved, the divergence of T^" [Eq. ([5])] is zero and the first law of 
thermodynamics is obtained projecting along f/^. We will do the same here but using 
qUIJij- as the projector. Although using Schutz's Lagrangian leads to work with a four- 
velocity normalized to —1 {qi = 0), the calculations in this section will be done with 
an arbitrary qi. The motivation is to obtain equations valid even when that condition 
is not satisfied. 

The results of this section will be applied to systems whose scales are much smaller 
than cosmological scales, so we take A^ = 0. Projecting T^.^ along c'^^U^ and using 

^On the other hand, the cosmological constant is important for the evolution of ip. 
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Eqs. (jSD, (fT7|) . and ([7]) together with Schutz's Lagrangian leads to0 



1 + 



^qi 



c 



dp 



(p + p) 



n 



dn= (p + p) + q2 -b^ dip . 



(18) 



It is convenient to express this relation in terms of the specific thermodynamic variables, 
v = — = —,u = pv = — , and = (V, U, and S are the total volume, total energy, 
and total entropy of a system containing particles): 



du 



+ pdv - v{p + p) + q2 -b^ dip + 



rfp= 0. 



(19) 



The first two terms are the only ones appearing in GR. The third term shows that ip 
formally plays the role of a new thermodynamic variable (changes in ip cause changes 
in the internal energy). Finally, the fourth term involves the pressure, which is not an 
independent variable (up to this point the independent variables are taken to be v and 
Ip). The expression is valid only for thermodynamic processes which do not involve 
heat transfer, while in a more general situation, with an amount dQ of heat being 
transferred, the equation becomes 



du + pdv - v{p + p) + q2 -b^ dip + 



dp = dQ . 



(20) 



This is the modified first law of thermodynamics in the VSL theory. 

Incorporating Caratheodory's principle it can be shown that there exists an inte- 
grating factor for dQ in Eq. fl20|) [25j . Defining ^ as the integrating factor and ds 
Eq. fl20l) becomes 

du + V 



dQ 
T ' 
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dp + pdv — v{p + p) + q2 — b^ dip = Tds 



(21) 



where s (identified with the specific entropy of the system) and T are two thermody- 
namic variables. To go on, it is convenient to introduce the function f{ip) = 1 + ^ = 

1 + ^e'^'^^. We take as the time of reference the present epoch, so that Uq = —1 (which 
is the usual normalization of the four- velocity) and Cq^ = c^l^^y- Then from Eq. (jS]) it 
follows that C = — Cq^ and hence f{ip) = 1 — e^^^. After choosing s, v, and ip as the 
independent variables and writing p in terms of them, an expression for the first law 
of thermodynamics involving only state variables is finally obtained: 



du = -{p + vfiiP)^] dv+fr- vf{iP)^^ ) ds + 



(22) 

(^{u + pv) + g2 - - '"fW^^ • 



^If Am 7^ but depends only on ip, the equation will be valid after adding the term —^^^j^d'i' 
to the right-hand side. 
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4.2 Equations of state 



The first partial derivatives of tlie specific internal energy can be obtained directly from 
Eq. fl22l) . There are three different equalities between the mixed partial derivatives 
which impose some restrictions on the functional dependence of p and T on ip: 

dv dT 

(23) 

dip dv^ 
dT dT 

9^ = "=^ + "."^. (25) 

with bi = ^ — q2 + b and 62 = y + q'2 — 

Several observations corresponding to different epochs in the history of the Universe 
show that the a variation has been very small (or zero) [2S], so the field ip must also 
be very small. Then, it makes sense to express the pressure as a power expansion in 

00 

p = J2pk{v,s)ip'. (26) 

k=0 

After replacing this series in Eq. and equaling terms with the same power oi ip a. 
recurrent formula for the coefficients pk is obtained: 

biv dpk{v,s) 

P..i(-,3) = ^^^, (27) 

where po{v,s) is the pressure as a function of v and s for = and therefore it is 
obtained from the usual theories in which c is constant. 

Working to first order in ip we arrive at the following expression for the functional 
dependence of p on f , s, and ip: 

p c:^ Po{v,s) + biv^ip . (28) 



5 Newtonian limit 

The Newtonian limit of this VSL theory can be obtained following the same steps 
as those used in GR. Although in the Newtonian limit time derivatives are negligible 
with respect to the spatial ones, special considerations are needed for the derivatives 
of ip- For example, if the spatial extension of a system is small compared to the scales 
associated with the spatial variations of ip (which are cosmological scales) and one 
is interested in following its time evolution, time derivatives will be more interesting 

■^The "ip dependence of the temperature can be treated in a similar fashion. 
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than the spatial ones. This will be the case in the following sections, where we will 
determine the evolution of the planetary radii and the luminosity of white dwarfs. We 
also remind that, as explained before, the cosmological constant is not included in our 
calculations. Finally, since the if) field must be very small, we will work to first order 
in ip. 

The condition of weak gravitational field allows to choose nearly Lorentz coordinates 
in which 

g^iv = ^l^lu + h^y, with < 1 . (29) 

Using Eqs. ^ and (jl]), the Lagrangian = p and the fact that aVoS/^ip is negligible 
(see the Appendix), the following gravitational potential equation is obtained: 

- V^/ioo = ^ Too , (30) 

where G = G {l + 2K+3a^ ) effective gravitational constant. In the nonrelativistic 

limit of GR /loo is identified with —2^ to arrive at Newton's gravitational potential 
equation, instead here we will go on working with /iqo due to the presence of the 
factor in that identification. 

The Euler's equations for this theory are obtained projecting Eq. ([5]) perpendic- 
ularly to with the projector ^f^" + U^U°' (here we have taken gi = from the 
beginning) : 

- ^Too/ioo,^ + P,^ + pU^'Ui^u + pUoU^goj,i = . (31) 
For quasistatical situations {Ui ~ 0) the equation becomes 

^TooV/ioo = Vp. (32) 



The hydrostatic equilibrium equation follows after combining Eqs. (l30l) and ([32 
For a system with spherical symmetry the result is 



dp G p{r)U{r) 



dr 



(33) 



where U{r) is the total energy inside the sphere of radius r and ^ must be understood 
as a spatial derivative at constant time. To obtain this equation it is necessary to 
consider c as a constant in the integral ^^J^ ^ dr. This can be done because the 
spatial variations of ip are negligible in noncosmological scales, as is demonstrated in 
the Appendix. 

6 Evolution of the radius of Mercury and a bound 
for the variation of c 

The radius of a planet is determined with the hydrostatic equilibrium equation to- 
gether with an equation of state and boundary conditions. The presence of ip in these 
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equations causes in general time variations of the radius. On the other hand, the actual 
change in size of several bodies of the Solar System have been estimated using different 
topographical observations. For Mercury there is a stringent bound: its radius has not 
changed more than 1 km in the last 3.9 x 10^ years ^27j. This fact will allow us to 
obtain a bound for the temporal variation of ip. 

In Sec. 6.1 we will show that the hydrostatic equilibrium equation is equivalent to 
another equation in which the temporal dependence resides only in the gravitational 
constant. Then it will be possible to use the results of the work of McElhinny et al. and 
give a bound for the quantity f{qi)ip, where f{qi) is a function of the parameters bi, q^, 
and g4 (the dot denotes derivative with respect to time). In Sec. 6.2 those parameters 
will be expressed in terms of the parameters originally defined in [Tn|. Finally, in 
Sec. 6.3 the bound obtained for f{qi)ip will be combined with a bound for the variation 
of the fine structure constant to obtain an upper limit for 



6.1 Transformation into a variable G theory 

To solve the hydrostatic equilibrium equation it is necessary to have an equation of 
state relating the pressure and density. It can be obtained replacing the corresponding 
equation of state for the constant c case in Eq. ( 1251) . For Mercury it is sufficient to 
work with a linear equation in p 



Po = Ksur(-^-l). (34) 

\ Po sur J 

The quantities with a "0" subscript correspond to ■?/' = 0, sur indicates evaluation in 
the surface of Mercury, and Ksur is the superficial compressibilitj0. Using Eqs. (1281) . 

2 

and the nonrelativistic expression po ~ — ^ (with mg the average mass of a 
particle), the equation of state to be used is obtained: 

p = K^ur ( — 1 ) - biKsuri^ . (35) 

\ Po sur / Po sur 

Replacing the equation of state (!35l) in Eq. (133!) . using the expression for ^ given 
in the Appendix and taking into account the definitions of q^ and q^ we get 

^ = -G e(^i+'?3+34-4)V' P»no(r)Mo(r) _ ^snr a ^b,^^ Pmo(r)Mo(r) 
dr Posur 2k + 3a2 ' 



where Pmo('") is the mass density for = and Mo(r) is the mass contained within 
radius r (also for ip = 0). The second term of the right-hand side is small compared to 
the first on^, so it can be multiplied by e^'^^~^'^^~^^'^ [which is equal to 1 + O^tp)] and 



In [27] a quotient between mass densities is used instead of , nevertheless in the Newtonian 
Hmit they are equal since Penergy ~ PmassC^- 

^Assuming &i,a = 0(1), the quotient between these terms is no larger than approximately . 
Since Ksur ~ lO^^ ^ and posur = PmOsur « 4 X 1021 [55], we see that ^ « IQ-^. 
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hence the equation becomes 



{r)Mo{r) 

dr ' 



with G the final effective gravitational constant 

a2 . K 



G = Gil + + 2h ^ . (38) 

\ 2k + Posur 2k + J 

The last term within the parenthesis depends on the particular properties of the body 
and could have, in principle, observational consequences. Nevertheless, its value is too 
small relative to the second (and body independent) term for its effect to be measurable. 

The method used to set a limit on the actual time variation of the radius of Mercury 
is based on the observation of surface features [27] . A homologous change (R oc M^/^) 
cannot be detected through this procedure since it scales all linear dimensions equally. 
Neither can be observed a variation i? oc r;, ( rf, is Bohr radius ), since a change in 
Tfe implies a change of all macroscopic dimensions in the same proportion. For these 
reasons it is necessary to make a change of variables: 

r , M(r) Mo(r) , , 

where M is the total mass of Mercury. Taking also into account that m oc c'^~^ and 
Tb oc c~'^, Eq. fl371) becomes 

^ = -G(t)^ p^^M*(r* 
dr* Tho r* 



_^(t)i:^ ^-°^'^^;;o^' \ (40) 



with 

Q(^l^ ^ Q g(6l+g3+<?4-4+2/3<7+2/3)i/.(i) j^^-j^^ 

and r^o the Bohr radius for if) = This equation shows that the VSL theory that is 
being studied is equivalent, with respect to the hydrostatic equilibrium of Mercury, to 
a theory in which the only "constant" that varies is G. 

Now we are in a position to use the results of [27]. The variation in the radius of 
Mercury (i?) produced by the variation in G can be parametrized as |^ 

1 dR 5 dG{t) 



R dxo G(t) dxQ 

or equivalent ly 

1 dR _ 5 dG{t) 
rH ~~ Git) dt 



(42) 



(43) 



where S is generally a function of G and M. Using models of Mercury, McElhinny et 
al. obtained the value 6 = 0.02 ± 0.005 for that planet. 



^We call "(5" the parameter called "a" in [57] to avoid confusion with the fine structure constant. 
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6.2 Specification of the parameters 

We will express the parameters bi,qi,q2,q3, and that have been introduced in this 
work in terms of the parameters q and b of the VSL theory. In our approximation 
M = NMn where N is the total number of nucleons. Equation ([7]) requires that 
N oc exp{—q2ip) while Eq. f[T^ implies that M oc exp{—bilj) and from the discussion 
in Sec. 2, we see that M oc exp [(g — 2 — q2)ip] and consequently q2 = b + q — 2. From 
Eq. ([8]) g3 = g — 2 + 3g = 4g — 2 while from Eq. ([9]) g4 = —q2 + q — 2. Also we have 
bi = qi/2 — q2 + b = —q + 2 and as we explained before we consider the gi = case. 
Finally from Eq. (14T|) we conclude that 



G{t) = G exp 



2 2 

&i + gs + g4 - 4 + -g + - ) ^ 



G exp 



— q — b ] ip 



(44) 



6.3 Bound for c/c 

Replacing the previous expression for G{t) in Eq. (H3|) one gets 

- i ^ ~ ± 5 X lO^^V"^ • 

with y~^ = 1/year. We have taken S = 0.02 ± 0.005 and ^ = ± 0.0004 [27], where 
AR corresponds to a time interval approximately equal to 3.5 x 10^ years. 

This result can be combined with bounds for a /a that have been obtained using 
atomic clocks. We can use e.g. the one obtained in Ref. 





(45) 



- = (4.2 ±6.9) X lO^^V^- 



(46) 



In this paper we will consider the 6 = case, which gives an upper bound for all 
non-negative values of b. Moreover, in the VSL theory - = qip, then Eq. fH5l) can be 



written as — ^V' = 0±5 x 10 ^^y ^ — Comparing this last equation with Eq. (jl6 
we see that — ~ ± 5 x lO^^^y^^ The conclusive result is that 



c 



0±2 X lO^^V^- 



(47) 



This can be rewritten as a bound for the adimensional quantity -ip' = Hq ^-ip after 
multiplying by the Hubble time {Hq~^) 0: 



V-' = ± 3 X 10 



-2 



(48) 



""We have taken Hq ^ ^ 1.5 x 10^" years. 
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7 White dwarfs luminosities and the scalar field 



White dwarfs are excellent objects to test any energy injection from a scalar field 
j30l [31]. This is due both to their low luminosity, as well as their extremely high 
heat conductivity, making all energy microscopically released to enhance the total 
luminosity. Most of them are adequately described by Newtonian physics and a zero 
temperature approximation, the latter hypothesis providing a polytrope type equation 
of state (EOS). 



7.1 The polytropic EOS 

Using the subscript to denote quantities without the presence of the scalar field, we 
write the polytrope equation as po = K^p^ . Given Eq. (128|) . we have 

P = Po{p) - Wlp—^, (49) 
P 

so we can write 

p = K^{1 - h^^)p\ (50) 

and we recover a polytrope equation of state with a new constant Kq ^ K = Kq{1 — 
biip'y). The Lane-Emden function with polytrope index (7 — 1)~^ still applies, and 
consequently the expressions for the radius, mass, and internal energy of the star will 
be the same in terms of the effective constant K p5j: 

3^GM^ (51) 

57-6 i? ' ^ ^ 

R = f'^'pt'^^'Ci, (52) 
M = 47rp^-^)/Vni^|^'(Ci)|, (53) 

where pc is the central density, (i is the first root of the Lane-Emden function, and we 
have defined / = 7^7/(477(7(7 — 1)). Unlike with the topographic analysis of planetary 
palaeoradii, the energy balance in the star is not invariant under the scaling (1391) . 
Thus, from Eq. (l37|l we see that the G for the internal energy of the star is G{t) = 
Gexp[{bi + qs + — 4)'0] = (5exp[(3g — b — 4)?/^]. On the other hand, the dependence 
of the radius R on ip can be obtained solving Eq. (l53l) for the central density and 
replacing it in Eq. ( l52l) [the ip dependences of the total mass and the compressibility 
were given in Sec. 6.2 and in Eq. fl5Q|) . respectively]. Finally Eq. flSTj) leads to the 
following dependence of the internal energy: 



E oc exp?/'/(g,6,7), (54) 

where 

,„,,„.3,-.-4.?l^l±M±i(i^. ,55) 
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To go on we assume that the star is in equihbrium in the sense that all the energy 
injected by the field ip is radiated away. Therefore the equation for the luminosity 
induced by the ip field (L^) becomes 

L^ = -E = -f{'y,q,b)Eij. (56) 



7.2 Comparison with observation 



We shall consider only white dwarfs well described by the nonrelativistic value 7 = 5/3, 
in which case f{q,b,'y) = y g — ^6 — ^. Unfortunately, only a handful of white dwarfs 
have well-measured masses, radii, and luminosities: indeed, these four or five stars 
are used to test the theory of white dwarfs, since the mass-radius relation requires 
exactly the same parameters we need to carry our comparison of theory and experiment. 
These stars and their properties have been reviewed in Ref. |32]. Table 1 shows the 
adopted values. We have excluded Sirius B from the sample, since relativistic effects 
are important in this case. 

To obtain the bound for ip we rewrite Eq. f l56|) as 



+ 



13 



a 



(57) 



E 136 + 14 ' 136+ 14 a 

Using again the upper bound P6|) for the present time variation of a and bounding 
by the observed luminosity of the white dwarfs (L) we obtain 



< — 



L 
E 



13b + 14 
3 



where 



and 



136 + 14 



^0 



En 



+ 



+ 



13 



a 



136+ 14 a 
13 



(5^ 



136 + 14 

(L/L^) 1 
{E/Eo)t^' 



3GMI 



1.1 X 10 



-14„ 



7 

is the would be internal energy of the Sun were it described by a Newtonian 7 
polytrope. Lq is the solar luminosity and 

To = 1^ ~ 1.32 X 10^?/, 



(59) 
(60) 

: 5/3 

(61) 



is the Kelvin-Helmholtz solar contraction time scale. 

The table shows upper bounds for ip in again for 6 = 0. 

Stein 2015B provides the strongest bound. Using a value for the Hubble time 
Hq^ oc 1.5 X 10^°?/ we obtain 



+ 2.1 X 10" 



(62) 
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Object 


M/Mq 


R/Rq 




E/Eq 




1^1 < 


Procyon B 
40 Eri B 
Stein 2015B 


0.602 
0.501 
0.66 


0.0123 
0.0136 
0.011 


5.8 X 10"^ 

0.014 
3.1 X 10"^ 


29.5 
18.5 
39.6 


1.5 X 10-^2 
5.7 X 10"" 
5.9 X 10^^^ 


3.3 X 10^13 
1.2 X IQ-^^ 

1.4 X 10-13 



Table 1: Data and bounds for selected white dwarfs (data from Ref. [32j ) . 



Comparing the expressions fj45l) and fl58|) we see that white dwarf physics provides 
the strongest constraints on the VSL theory near the present epoch for almost all 
values of the b parameter, except for those near h = —14/13 which make Eq. fl58|) 
uninformative. It is also clear that combining both bounds fllS]) and flSSl) we can 
obtain a bound for independent from the value of b (although less strong than the 
bound given for 6 > 0): < 2.2 x IQ-^^y^^. 

8 Conclusions 

We have obtained the equations that describe a perfect fluid in the nonrelativistic 
limit and the first law of thermodynamics in the context of the covariant VSL theory 
proposed by J. Magueijo. We showed that the field ip can formally be considered as a 
new thermodynamic variable and we also showed how to obtain the equations of state 
in the VSL theory when the corresponding equations for constant c are given. 

The nonrelativistic hydrostatic equilibrium equation has the usual form with the 
gravitational constant G replaced by an effective constant. The different variables 
(pressure, mass, mass density) depend on the ip field and so the radius of a planet 
should vary in time. Using bounds for the variation of the radius of Mercury and the 
fine structure constant we have set limits on c/c: c/c = ip = ±2 x 10~^'^y~^ (valid 
for positive values of the b parameter). The most interesting thing about this result is 
that it gives a bound for ip, whereas the known limits for the variation of a and e lead 
to bounds for the product q'lp- 

Under the same Newtonian approximation we obtained the dependence of the lu- 
minosity of a white dwarf on the time variation of the scalar field. The bound obtained 
is more stringent than the planetary radius bound by an order of magnitude. 

The 6 = assumption suggests a null coupling of the scalar field with matter. 
However, the q assumption implies a quantum coupling between ip and matter, 
not explicitly shown in the action ([2]). Of course this and other issues such as the 
microscopic origin of the energy exchange between the scalar field and ordinary matter 
as well as whether all the energy injected by on a star is radiated away or not, deserve 
further work. This we leave for future communications. 
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9 Appendix: spatial and temporal behavior of ijj 



After some simplifications valid in the Newtonian limit, the equation for ip (in the case 
a 7^ 0) becomes: 

= ^""^ aT . (63) 
^ c4(2K + 3a2) ^ ' 

Since is small, can be replaced by Cq in the right-hand side of the equation (this 
is the first step of an iterative process of resolution): 

= 4,/""f o 2, «^ • (64) 

T can be written as the sum of a spatial part (corresponding to Mercury) and a 
temporal one (cosmological term). Then ip can also be separated into spatial (ips) 
and temporal (ipt) components. The metrics to be used are the quasi- Minkowskian and 
the cosmological (Robertson- Walker) for the spatial and temporal parts, respectively. 



9.1 Spatial behavior 

Using that T ^ p = energy density of Mercury, the spatial part of Eq. 0641) can be 
written: 

^"^'-r^drV dr )~ c^f 2«: + 3a2^- ^^^^ 



(66) 



After integrating one arrives at 

dips _2G a U{r) 
dr Cq 2k + Za? 



This is one of the formulas that has been used in our work. It can help us to obtain 
an estimation of the variation of if) inside Mercury: 

2a 4 1 / G A 

where pm and R are the average density and the radius of Mercury, respectively. The 
first factors are presumably 0(1) so 



G 



Aijj ^ —p-^W ^ IQ-'' . (68) 
Cq 

This number is small compared to the bound that is obtained for the variation of ijj 
in the last 3.9 x 10^ years (the relevant period of time for our study of Mercury) and 
besides the temporal and spatial behavior of this field can be separated. This justifies 
the steps of our analysis in which ip was considered constant inside Mercury. 
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9.2 Temporal behavior 

The temporal part of Eq. (IMl) is 



where Pmi7 is the average density of the Universe. Using the Robertson- Walker metric 
it becomes 

„ „n , SttG a hit) ■ 

where a(t) is the scale factor of the metric and a is one of the parameters of the 
VSL theory. In the way towards obtaining the gravitational potential equation in the 
Newtonian limit, the equation = — |i? = Tq° + a (Vq V^?/' — O^) appears. 
We want to show that VoV'^ip is negligible compared to the order ip term in the Taylor's 
expansion of Tq°. To do this it will be demonstrated that each term appearing in 
Eq. (170|) can be neglected: 

• First term of Eq. fITOl) : This term is negligible since Tq° ^ PMercuryC^ and ^""^ ~ 
10-30. 



PMc 



Second term of Eq. fITU]) : 

a(t)^^^AV.^ 1^ (71) 

where Aip represents the change of in a time interval At. Taking At ^ Ax 10^ 
years (this is the time interval for which we have a bound for the variation of the 
radius of Mercury) , we get 

3^ij ^ 10-^^ Ai: cm-^ (72) 
a[t) 

This quantity must be compared with 

(-4) pij ^ 10"^ X — X p^x ij ^ lO-^V cm-2 , (73) 
and hence we see that the second term of Eq. (ITOll is also negligible. 
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